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Several puzzling regularities concerning the low temperature excitations of glasses are quantita- 
tively explained by quantizing domain wall motions of the random first order glass transition theory. 
The density of excitations agrees with experiment and scales with the size of a dynamically coherent 
region at T g , being about 200 molecules. The phonon coupling depends on the Lindemann ratio for 
vitrification yielding the observed universal relation l/X ~ 150 between phonon wavelength A and 
mean free path I. Multilevel behavior is predicted to occur in the temperature range of the thermal 
conductivity plateau. 

PACS Numbers: 66.35. +a, 64.70.Pf, 66.60.+f 
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Decades ago, measurements of the heat capacity and 
thermal conductivity of glasses at cryogenic temperatures 
revealed the presence of excitable degrees of freedom not 
present in perfect crystals 0]. These could be described 
as two level tunneling systems whose energies and tun- 
neling matrix elements were randomly distributed [pH. 
Coupling the tunneling systems to phonons explained the 
thermal measurements and also predicted novel physical 
effects, such as the nonlinear absorption of sound and a 
phonon echo, which were later observed |HS- 

The nature of the tunneling entities has remained ob- 
scure. Thoughtful experimentalists and theorists have 
noticed puzzles when the model quantitatively confronts 
experimental data For example, the entropy con- 

tained in these excitations is much less than the residual 
entropy frozen in at the glass transition. Yet, the density 
of two level systems varies only modestly from material 
to material. There is also a mysterious nearly universal 
relation between the density of the two level systems and 
their coupling to phonons which can be deduced from the 
observation that the mean free path of phonons is about 
150 times their wave length at low temperature 0. If 
the two level systems arise from the motions of highly 
localized specific configurations of atoms, as in impurity 
doped crystals, instead of such a universal relation we 
would expect significant variation with the glass's chem- 
ical composition. These facts lead Yu and Leggett ||, 
as well as others |^,|l0) , to investigate the possibility that 
the experimentally observed excitations are really highly 
renormalized collective excitations of a system of micro- 
scopic tunneling entities that interact strongly through 
the exchange of phonons. While such a coupling seems to 
be present, manifesting itself in spectral diffusion of the 
two level entities 11 , quantitative calculations based on 



the interacting model suggest that thermodynamic man- 
ifestations of the interaction should be confined to ultra 
low temperatures |H],[l2| ■ This scenario then has not yet 
explained the observations which called it forth. 

Here we explore an alternative view of the quantum ex- 



citations of a glass. Rather than regarding the tunneling 
entities as extrinsic we quantize the excitations that are 
responsible for the activated dynamical events in a liquid 
which slow as the glass transition is approached. These 
excitations are mostly frozen in at the liquid glass tran- 
sition. Many aspects of these activated motions can be 
understood using the random first order transition the- 
ory of glasses |0| . This theory starts from some exactly 
solvable mean field glass models showing one step replica 
symmetry breaking p^ , p^| . The picture of activated mo- 
tions goes beyond mean field theory by considering en- 
tropic droplets The theory quantitatively explains 
for a wide range of substances both the barrier heights 
p7[ and nonexponentiality of relaxations observed 
near the glass transition. In this picture, a viscous liquid 
or glass consists of a mosaic of frustrated domain walls 
separating regions of energetically less frustrated mate- 
rial. Each mosaic cell resembles a local minimum of the 
free energy, an "inherent structure" if you will fl9|| . The 
typical size of the cells is nearly universal and does not 
vary much with the composition of the glass because it 
depends on the near universal value of the maximum vi- 
brational amplitude sustainable by the glass - the Linde- 
mann ratio. The mosaic length scale depends not on the 
molecular character of the glass but on its preparation 
time scale, logarithmically. 

Large scale motions such as these domain wall move- 
ments have usually been discounted as a possible origin 
of the two level systems because tunneling amplitudes 
decrease with the number of independently moving en- 
tities and with the distance over which they move (see 
ref. p0}-|2^] for exceptions). We will show that occasion- 
ally in the amorphous solid such collective tunneling can 
occur with a significant amplitude owing to the extraor- 
dinary multiplicity of tunneling paths available for such 
an entity and to the possibility of achieving a sharp local 
resonance. 

The density of excitations available through collective 
tunneling can be calculated using an argument based on 
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the classical density of states of systems with one step 
broken replica symmetry p^ |. The result depends only 
on the glass transition temperature T g and the mosaic 
size, and explains the magnitude of the heat capacity. 
The coupling of these tunneling excitations to phonons 
is set by T g itself, again because of the universal Linde- 
mann ratio. For resonant scattering, the material depen- 
dent factors cancel to yield the observed relation between 
phonon wavelength and mean free path. Our arguments 
show that tunneling will appear to be two state like over 
the wide-range of temperatures where echo experiments 
have been performed, but that the two state character 
breaks down at higher temperatures where a plateau in 
the thermal conductivity is observed H ] and where sin- 
gle molecule experiments already give evidence for devi- 
ations from two-stateness |25| . 

First, a simple argument for the density of excita- 
tions to set the stage. The motions above the glass 
transition temperature T g are rearrangements of finite 
sized cooperative regions from one local free energy min- 
imum to another ||l7). The free energy cost to create 
a new minimum as a function of the droplet radius is 
F(r) = 4ira(r)r 2 - ^-r 3 TS c . Here, a{r) = a^a 1 / 2 ^' 2 
is a radius dependent surface tension whose form fol- 
lows from a renormalization group calculation and is 
caused by the wetting of the interface between two par- 
ticular low energy configurations by other possible ar- 
rangements, do is surface tension at the molecular length 
scale a. s c is the configurational entropy per unit vol- 
ume, which favors creating other minima. In terms of 
the number of particles in the droplet N = 4^(r/a) 3 , 
we have F(N) = j^/N — Ts c N, where 7 = 2v / 37rcroa 2 , 
and s c (T) is configurational entropy per particle. An ap- 
proximate density functional calculation jl7| gives 7 = 
I \/3tt k bT ln(a£, a 2 /ire), where a_La 2 ~ 10 2 , the inverse 
square of the Lindemann ratio, hardly varies from sub- 
stance to substance. The maximum of the free energy 
F* = j 2 /4Ts c is reached at N = (^/2Ts c ) 2 giving the 
typical motional barrier. T g is set by the quenching time 
for the glass r = Toe pt l kBT a ; where To is a molecular 
time scale. At T g the system breaks up into a mosaic 
of regions of size £, where F(N*) — 0, giving N* ~ 190 
and £ = N* 1 ^ 3 a. £ only logarithmically depends on r 
and is about ~ 5.7 molecular radii, for quenching times 
of hours, independent of molecular composition. 

Each mosaic cell resembles a finite size mean field sys- 
tem at Tk i.e. where the entropy vanishes. The den- 
sity of minima n(e) for any system experiencing an en- 
tropy crisis is: n(e)de = k c Hp e e ' kBTg . The proportion- 
ality constant l/ksTg guarantees that there is only one 
absolute ground state: J_ den(e) — 1. The total den- 
sity of states in the sample per unit volume is therefore 
N(e) ~ k J r ^ a e e / fcgTg . We do not distinguish between 
the energy and the free energy of basins because the 
zero point energy and vibrational entropy vary little from 



basin to basin. At low energies, one thus gets for the con- 
ventional density of excitations P = ,, A-gg- . For silica 

with T g ~ 1500 K and a ~ 3.5 A, this gives ~ 6 • 10 45 J" 1 
m -3 , a value typical for many glasses pi) . Apart from 
the T g variation, which has been noted experimentally 
(see e.g. [^6|), the density is seen to be nearly universal 
for glasses made with the same quenching rate. By inte- 
grating the flat distribution P up to T g one finds that the 
low level excitations accessible at cryogenic temperatures 
would account for less than (a/£) 3 ~ 1% of the residual 
entropy at T g . 




FIG. 1. A schematic of a tunneling center is shown. £ is 
its typical size, is a typical displacement of the order of 
the Lindemann distance. Red contour illustrates a transition 
state size. 



A more complete argument demonstrates that reso- 
nant tunneling consisting of small displacements of atoms 
within a region of size on the order £, as illustrated 
in Fig.^, can indeed occur and gives the same excita- 
tion density as the simpler argument. We first exam- 
ine the existence of resonant levels. Consider a region 
from the system quenched below T g which forms a par- 
ticular local minimum. An excitation corresponds to in- 
troducing a different minimum structure, encompassing 
N molecules. The internal energy will be on average 
Acp(T ff — Tk)N ~ T g s c N, where Ac p is the configura- 
tional heat capacity. On top of the internal energy, there 
is a contribution from an interface energy 7V^V due to 
the mismatch with the configuration of the surrounding 
material. The distribution of excitation energies of those 
"droplet" configurations should be gaussian p5[ ], giving 
for the number of states higher in energy by E than the 
original one n N (E) ~ exp {s c N - H^W^lI }, 

SE 2 is the variance in energy per particle at the glass 
transition temperature. s c is the residual entropy that 
would be frozen at T g [|l7| . The entropy crisis for the bulk 
fixes the glass transition temperature by 8E 2 = 2s c T K . 
The number of states at small excitation energies grows 
with the excited region's size. Yet below a critical size 
the domain wall energy prevents resonance. The number 
of nearly resonant levels only becomes on the order one 
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for N — ("f/TgSc) 2 N = N*, as before. Again, excita- 
tion energies are distributed according to an exponential 
e E/k B T g ^ gi v i n g the same P. The multiplicity of nearly 
isoenergetic i.e. resonant configurations intrinsically fol- 
lows from the non-equilibrium freezing at T g . 

The enormous multiplicity of states at T g not only al- 
lows resonance but also yields myriads of possible tunnel- 
ing paths from the initial configuration to the resonant 
one with a droplet of size N* embedded in it. Each path 
is a connected sequence of droplet-like configurations in- 
volving rearrangements of a region containing a growing 
number TV < N* particles. The tunneling amplitude will 
be a sum over all these paths, much like the partition 
sum for a random directed polymer |p7fl with the weight 
of each path being the exponential of its action in units of 
h. If h is large, the path sum will be dominated by a large 
number of paths. These would go through a nearly con- 
tinuum of paths as shown in Fig.|^. This situation would 
correspond with a "quantum melted glass" Q . On the 
other hand for small h, the smallest action path only will 
contribute giving a tunneling clement e ~ s " li ™/' i . This ac- 
tion varies from one resonant pair to another and will be 
exponentially distributed over a range scaling with ksT g , 
giving a power law distribution for the tunneling ampli- 
tude. Finding the precise statistics of the lowest action 
path from the direct polymer analogy is quite complex 
since the tunneling distribution depends on detailed (cor- 
related) statistics of the energy variations all along the 
tunneling paths. We obtain a sensible approximation by 
recognizing that the action will be crudely proportional 
to the largest barrier encountered as the virtual tunneling 
droplet grows, as would be true for an inverted oscilla- 
tor potential. To find the statistics of this barrier at any 
value of N we must find the distribution of the lowest 
energy droplet state at N (notice we neglect correlations 
here). This distribution is of the same form as the one 
in Qn{E), except the variance is twice larger now be- 
cause the energies of both initial and the transition states 
may fluctuate since the system can choose where to be- 
gin to tunnel: fl N (V) ~ exp {s c N - [V ' {T ^^ )? }■ 
Generally, where a resonant state exists the tunneling 
path will start by rearranging high energy local con- 
figurations into ones with internal energies near Ex- 
The lowest barrier likely to be encountered occurs when 
^Af(Knp) — 1) as usual in the extreme value statistics, 
giving Vmp — "/VN — (2\/2 — l)T g s c N . V as a function 
of N is shown as solid black line in Fig.^|. The maximum 
value of V mp is proportional to the activation barrier at 
T g , F* but smaller due to the fluctuation of the initial 
energy: V max = F^/{2^2 - 1) ~ 26T g s c . The maxi- 
mum occurs at N* = N Q /(2y/2-l) 2 ~ 14. This is small 
enough that we can expect system dependent corrections. 
There will generally be states below V mp which can allow 
tunneling as in the green line shown in Figj^. The dis- 
tribution of barriers (and therefore actions) below V rnax 



follows from f2jv(V) giving an exponential distribution 

proportional to 0(y*) ~ exp j — 18 ■ s c + j. The 

chances to be able to tunnel to a state of precisely the 
minimum resonant size is suppressed by a factor e~ 18Sc . 
But to find a state which is simultaneously resonant and 
within easy tunneling requires encompassing a region 
with only an 18 additional molecules, less than a single 
layer. Hence, any region of size ~ 200 molecules will have 
a nearly resonant tunneling state within range T g . Tun- 
neling involves simultaneous motion of all the atoms in 
the droplet, and might have a high effective mass or even 
be damped owing due to the complex rearrangements in- 
volved. It is hard to rule out the latter possibility, but the 
actual effective mass is low since moving a domain wall 
over a molecular distance a in an (imaginary) tunneling 
time r involves displacing individual atoms only a Linde- 
mann length d^. The kinetic energy associated with this 
motion is M w (a/r) 2 — N w m(dL/T) 2 , where N w ~ (£/a) 2 
is the number of molecules in the wall and m is the 
molecular mass. Thus the mass of the wall M w is only 
m{^/a) 2 {d L /a) 2 . Using (f/a) ~ 5.8 and {d L /a) 2 ~ 0.01 
gives M w ~ m/3. Computational studies of similar 
multiparticle tunneling events |2p| , ^l|| are consistent with 
the low mass obtained here. As we shall see below, 
(dhja) 2 ~ ksTg/pc 2 , where p is mass density and c s 
is the speed of sound. This gives M w ~ (£/ 'afhsTgl 'c 2 . 
It follows that the frequency of motion at the barrier 
top d* = — d 2 V/dr 2 /M w ~ 1.6(a/£)w£>, expressing V as 
function of droplet's radius r = a(3A r /47r) 1 / 3 . 
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FIG. 2. The black solid line shows the barrier along the 
most probable path. Thick horizontal lines at low energies 
and the shaded area at energies above the threshold repre- 
sent energy levels available at size N. The red and purple line 
demonstrate generic paths, green line shows the actual (lowest 
barrier) path, which whould be followed if %ur < ksT /2n. 

Since the tunneling matrix element A for a path with 
barrier is proportional to e ^ 7rVt / h ^ t we obtain a 
rather flat distribution for the tunneling exponent log A, 
as is generally used to fit experiments. 

We now determine the size of the coupling of the do- 
main motions to the elastic (single polarization) strain 
field V0, which has potential energy density pc 2 (W(f>) 2 /2. 
At low temperature the phonon wavelengths are long, so 
one can describe the interaction with a standard point 



3 



like term a\ gV0, where <j z is an operator of the tunnel- 
ing degree of freedom. The tunneling entities first come 
into existence at a temperature Ta, somewhat higher 
than T g , where mechanical stability of local minima to 
thermal vibrations is achieved [ ft3[ . The phonon ener- 
gies at the microscopic scale and their coupling to the 
defects will be comparable and on the order of T g , giving 
pc 2 a {(V(j)) 2 )a 3 ~ (gV0 a) ~ k B T g . This sets th e cou- 
pling g at the molecular scale g ~ \J pc 2 a 3 ksTg. The 
coupling to the extended defects is weakly dependent on 
their size. To see this, in the continuum limit, we sep- 
arate the total elastic deformation tensor iiy into con- 
tributions <f>ij and ±dij/2 due to phonon displacement 
and the tunneling motion respectively. The difference in 
energy between the defect configurations in the presence 
of a (longitudinal) phonon is then pc 2 <pu J d 3 rdu, where 
the integration covers the droplet. The coupling is there- 
fore proportional to a surface integral J dsd(r), where 
d(r) are the tunneling displacements at the edge of the 
droplet. These are random and of order (a/£)dz, because 
the inelastic displacements decrease from ~ dj, in the 
center of the droplet to zero outside. The integral is of 

order a 2 V N* 2 ^ 3 (a/^)dL, therefore g ~ pc 2 s a 3 dL- Using 
dhja ~ V</> at T g , one still gets g ~ yj pc 2 a 3 ksT g . 

The resonant scattering from the tunneling of mosaic 
cells gives l^j p (w) = 7r^-wtanh ( 2 ^ T ^ as the inverse 
mean free path of a phonon with frequency to p9| , p4[ . 
Thus, Xds/lmfp - 3Pg 2 /pc 2 . Combining P ~ l/k B T g e 
with the expression for the coupling constant, one obtains 
Imfp/^dB — {£,/a,) 3 ~ 10 2 . This ratio depends only on 
£/a, independent of molecular composition. It is a geo- 
metrical factor reflecting the relatively low concentration 
of cooperative regions in a supercooled liquid frozen on 
quenching, an almost universal number within the ran- 
dom first order glass transition theory Q . 

At high temperatures the domain wall motion will be- 
come noticeably multilevel. Ignoring damping, at a tem- 
perature T' ~ hco* /2-Kks — (o/^)Td/27t, the wall motion 
will typically be classical. This temperature lies within 
the plateau in thermal conductivity [?|. Damping, which 
becomes considerable also at these temperatures, should 
lower this estimate, as also will the fluctuations in barrier 
height. A multilevel system will more effectively scatter 
phonons, which would cause the plateau. Consistent with 
this, single molecule studies of spectral diffusion of dyes 
in polymer glasses at these temperatures reveal spectral 
trails that wander fl25fl, as expected for domain walls in 
crystalline materials ]30|. 
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